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Abstract 

In this paper we prove the existence of phase transitions at finite temperature for 0(n) 
classical ferromagnetic spin models on infrared finite graphs. Infrared finite graphs are infinite 
graphs with lim m ^ + Tr(L + m) _1 < oo, where L is the Laplacian operator of the graph. 
The ferromagnetic couplings are only requested to be uniformly bounded by two positive 
constants. The proof, inspired by the classical result of Frohlich, Simon and Spencer on 
lattices, is given through a rigorous bound on the average magnetization. The result holds 
for n > 1 and it includes as a particular case the Ising model. 
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1 Introduction 



The study of statistical models on graphs requires the introduction of new techniques and 
concepts with respect to the well known case of lattices. The lack of translational invariance 
in general implies inhomogeneity in local variables and makes useless the introduction of such 
a powerful tool as Fourier transforms in space. In addition, due to the absence of a natural 
definition of dimension, the statement itself of general results and theorems is rather problematic. 
The key point is to find a connection between geometry and physical properties. This problem 
can be studied in the framework of algebraic graph theory. However, phase transitions only occur 
on infinite graphs and the algebraic theory of infinite graphs is a very recent field of research 
in mathematics Up to now only few classical results on lattices have been extended to 

graphs. 

In this paper we deal with the generalization of the proof of existence of spontaneous mag- 
netization at finite temperature for ferromagnetic models on lattices in d > 3. This result was 
proven in 1976 by Frohlich, Simon and Spencer || by a general approach based on Gaussian 
domination of correlation functions in the infrared regime. Such approach, as well as the for- 
mulation of the result itself, is deeply related to the translation invariance of lattices. Therefore 
here we modify the mathematical techniques and we use an alternative definition of the order 
parameter. In particular the concept of dimension is replaced by the asymptotic behavior of the 
Laplacian spectral density at low eigenvalues, a choice which is meaningful also in many other 
contexts. The structure of the paper is as follows: in section 2 we recall the basic definitions and 
theorems of graph theory to be used in the proof; in section 3 we give the fundamental known 
results about statistical models on graphs and state our theorem on the magnetization bound; 
finally in section 4 we present the proof. 

2 Some mathematical properties of graphs 

Definition 1 A graph G consists of a countable set of vertices (points, sites) V(G) = {i,j, k, . . .} 
and of a set E(G) of unordered pairs of vertices. The generic element of E(G) is called a 
link and its vertices are said to be nearest neighbors (adjacent). IfV{G) is finite, G is called a 
finite graph and we will denote by N the number of vertices of G. 

Definition 2 A path connecting two vertices i and j is an alternating sequence of vertices and 
links [i,(i,k),k,(k,h),h... ,w,(w,j),j]. 

Definition 3 A graph G is connected if, given any two vertices i,j £ V(G), it exists a path 
between i and j. 

Here we will deal only with connected graphs. The topology of G is algebraically described by 
its adjacency matrix. 
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Definition 4 The adjacency matrix ofG, Aij, is given by: 

r i if (,,) , E(G) 

{ otherwise 

Most physical models on G can be defined in terms of the Laplacian matrix Lij: 
Definition 5 The Laplacian matrix of G, Lij is given b: 

T J ij — %iO~ij -^-ij (2) 

where z% = J2jev(G) Aij is called the coordination number of i. 

In the following we will consider graphs with coordination numbers bounded from above: z% < 
Zmax G N Vi G y{G). If G is infinite, i.e. V(G) is infinite, can be considered as the 
representation of the Laplacian operator L acting on £2 (V (£?))■ Under the previous conditions 
for G, it can be shown that L is symmetric, non negative and bounded | |10| |. A relevant property 
of L we will exploit in the following is the Schwinger-Dyson identity J?]]. Let us consider a 
diagonal bounded operator rjij = rjiSij with rji GC such that (L + ry)~- exists. Then: 

J2 th(L + 17)^ = 1. (3) 

ieV(G) 

This identity follows from the property J2i Lij = 0. 

The graph G is naturally provided with an intrinsic metric defined by the chemical distance 
on G. 

Definition 6 The chemical distance r^j G N between two vertices i and j is the number of links 
in the shortest path connecting i and j. 

The intrinsic metric is the fundamental tool to define the thermodynamics! limit on a infinite 
graph G, allowing us to introduce the generalized Van Hove spheres. 

Definition 7 The Van Hove sphere Si jr of center i and radius r is the subset ofV{G), Si^ r = 
{j G V{G)) \rij < r}. We will define N^j. as the number of vertices in Si^ r . 

The behaviour of iVj r as a function of r characterizes the growth of G. In particular if 

Air c < Ni r < Bir c where c, Ai and Bi are positive constants Vi. (4) 

G is said to have polynomial growth and we define the connectivity exponent c to be the inf of 
the set of all c satisfying (||). Here we will consider graphs with polynomial growth. Indeed this 
condition is sufficient to guarantee the convergence of the restriction Oj /r to Si >r of a bounded 
operator O defined on G to O itself for r — > 00. 
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Definition 8 Let G be an infinite graph satisfying all previous conditions. Given a complex 
function of the vertices 4> : V{G) — > C we define the thermodynamic limit average <f> of <j> as: 



to* aT" £ (5) 



(/> can 6e easily shown to be independent of i. 

In this work we will use also on finite graphs the symbol 4> to indicate the average over all 
vertices. We will define the average on a subset V\ C V(G): 

Xv^^l™^- Xv x (i)<l>i- (6) 
where Xv x {^) ls the characteristic function of V\. We define the measure fJ,(V\) of V\: 

KVx)=m- (7) 

More generally, the thermodynamic limit average of a k variable function $ : V((G)) k —>C is 
defined by: 

1= run _L v *i lt ..,i h . (8) 

' ii,...,ifc£Ai,r 

We define the thermodynamic limit average trace of an infinite matrix B^ as: 

Tr5 = b (9) 
where hi = Ba, and its restriction to a subset V\ C V(G): 



Try, 5 = xv A t (10) 

The large scale topology affecting the critical behavior of statistical models on graphs can 
be characterized by the properties of the trace of (L + M) _1 , where M« = m 8™ is a real and 
positive diagonal matrix. 

Definition 9 An infinite graph is infrared finite if 

lim Tr(L + M)" 1 < oo. (11) 
An infinite graph is infrared infinite if: 

lim Tr(L + M)" 1 = oo. (12) 

m— >0+ 

It can be shown that if we substitute the constant matrix M with M-j = m hi 5ij where 
< e < bi < K the infrared behavior for m —>■ + does not change ||]. 

In dealing with statistical models the adjacency matrix Aij is often generalized to the fer- 
romagnetic interaction matrix Jij: 
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Definition 10 

j j j Jij with < e < Jij < J if = 1 

if^ = 1 j 

In i/iis case the Laplacian generalizes to: 
where Z{ = J2j Jij ■ 

The generalized Laplacian has the same properties of the Laplacian associated to A: it is 
symmetric, positive and bounded, it satisfies the Schwinger-Dyson identity and finally it has 
the same infrared behaviour iBjl. 



3 O(n) classical spin models on graphs 

An important class of classical statistical spin models on graphs is defined by the Hamiltonian: 

H ' = -\ E Jij#i'fj-H- E *i ( 15 ) 

i,j£V(G) ieV(G) 

where Jij is a ferromagnetic interaction matrix defined on the graph G satisfying conditions 
fli~3| ) and are n-dimensional real unit vectors = (aj, ... , of) defined on each vertex i and 
satisfying the constraints: 

1?1 = 1 Vi. (16) 

For n = 1 H' defines the Ising model which is invariant under the discrete symmetry group 
Z2, while for n > 2 H' represents an 0(n) model with continuous symmetry. Finally h = 
(h, 0, . . . , 0) is an external magnetic field coupled to a^j. In the following we will set h > 0. 
Due to (111) H' is equivalent up to an additive constant to: 



H = ^E^(^-^) 2 -^-E^ = ^E L ^-^-^-E^ ( 17 ) 

ij i ij i 

where the Laplacian operator defined in ( |i~4| ) has been introduced. 
The Boltzmann measure up h(~<?) is given by: 

d/i Aft (^) = n^e _/3H (18) 



where (3 = 1/fce T, with ~&i satisfying the constraints (|lq). The order parameter of the model 
is the average magnetization in the h direction: 



dnp jh {^) a 1 

M(J3,h) = *— r . (19) 

dfJ>p,h(~rf) 
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On lattices, in the thermodynamical limit for n > 2 it is known that the limit for h — ► of the 
order parameter vanishes for T > in the infrared infinite case (i.e. if the Euclidean dimension 
d is 2 or 1) In the infrared finite case (d > 3) for low enough temperatures the limit is 
positive P], implying the existence of phase transitions. On graphs we have lim/i^o M(/3, h) = 
for every j3, if for each positive measure subset V\ C V(G) we have 

lim TtyJL + b)- 1 = oo . (20) 

For Ising models on lattices the limit of the order parameter vanishes for all temperatures only 
in one dimension. There are no analogous results for generic graphs. 

In this work we will show that, for the classical spin models ( |T5| ) on infrared finite graphs, 
below a certain temperature M((3, h) satisfies the bound: 

lim M(p, h) > k > 0. (21) 

h— »0 

This result for 0(n) models is not the exact inversion of Q. Indeed one should prove that if 
there is a positive measure subset V\, for which the limit in ( p0|) is finite, then the bound ( ^l|) is 
satisfied. Here the proof is given only for the case where the limit in ( p0[ ) is finite for V\ = G. A 
complete inversion of j| will be given in a following paper Q. We remark that this is the first 
result for the Ising model on a generic graph. 



4 The magnetization bound 



In this section we will prove the magnetization bound ( |2l| ) for the classical spin model (15) 
according to the following strategy: 

(A) first of all we will prove some useful inequalities for the Laplacian on finite graphs; 

(B) then we will introduce for the constraint (^) an integral representation allowing us to 
perform a Gaussian integral with respect to the spin variables "a^; 

(C) we will take the thermodynamic limit of the inequalities (A) and in the expressions for 
M and for the global constraint; 

(D) exploiting a saddle point technique for large f3, we will apply the inequalities in (C) to 
M getting the bound (|2T|), for infrared finite graphs. 

0.3truecm As for the point (A): 

Lemma 1 The following inequalities hold for finite graphs, with h, L j = h dij and a = a% <5jj with 
Oi{ G R Vi: 

1 > Re (|E(^ + /i + ^ 1 ) >te(^ y £(L + h + ia)J*\. (22) 



°~nT, a ^ L + h ~ 2 «) J fc 1 ( L + h + ia )kj a i ^ Re E a i( L + h + ia )ij la i ] ( 23 ) 

ikj 
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< ^ReTr Vx (L + h + ia)" 1 ) < -^ReTry, ((L + h)' 1 ) < ^ReTr ((L + h)' 1 ) (24) 

^\JmTr Vx (L + h + ia)- l \ < ^Tr y , (L + h)" 1 < ^{L + hy 1 (25) 
where V\ is a generic subset V{G) and 

Tr Vx B = Y,B iiX v x (i) , (26) 

i 

( \ V2 

^ £«? < /(a) - J] + h- ia)- k \L + h + iafcfe (27) 

i \ ikj J 

where f{a) < l / N ^^af (l max + h) 2 + af) and l ma x is the maximum eigenvalue of L. 

The proof of these inequalities will be given at the end of this section. □ 
As for (B), we prove the following: 

Lemma 2 On a finite graph the expression of the average magnetization ( |7^ j and of the con- 
straint: 

f T^d^lf) 
1 = 

J dfj,f3 )h (~d>) 

can be written in the following way: 

M(P, h) = ^J d mh (a) A J2(L + h + ia)~ k j (29) 

kj 

d^,h( a ) = II date* 3 ** (31) 



(28) 



where: 



iS^h = --Tr (ln(L + h + ia)) + | I i ^ a< + /i 2 ^(L + ft + ia)^. 1 



and 



Z = J dfip, h (ai) 

Let us write the constraints (|I^) using the complex integral representation of the delta function: 

5(x 2 ) = / da exp(— iax 2 /2 — ex 2 /2) 



where e is a real arbitrary constant. We will put e = h(3. Substituting the expression for 5 in 
(ID and (||) we obtain: 



M(P,h) 



J dlf da ^exp {-P{\Y,Uf^i ■ Hi ~h^}) ~ E(^f(^ 2 - 1) - 
J dlf daexp ( -0(1 J2 ■ *i -hj^aj) - £(i |(^ 2 - 1) - M^ t ?) j 

\ ij i j / 



and 

/ ^ da ^exp ( £ • *i " * E **) " f(^ 2 " 1) " 

ij i i 

| d7? daexp ( • ^ -^^) - E(*f(^ 2 " 1) ~ ^)) 

\ *j i « / 

where da = Ylidai and dl? = Ylidlfi. Substituting in both integral a\ with (3ai and with 
CT^i/y^, we can perform the Gaussian integral on the variables o*j, obtaining in this way ( |29| ) 
and ©. □ 

Now we consider an infinite graph G and using the thermodynamic limit, we will extend to 
G the Lemmas (1) and (2). This is point (C) in our plan. 



Lemma 3 For an infinite graph the inequalities (\28i), ft23l), j$4Q, (\2di), t\2j) become: 



1 > Re [h(L + h + ia)" 1 ] > Re [h\L + h + ia)~ 2 j . (32) 

< a(L + h- ia)- l (L + h + ia)- 1 a < Re (a(L + h + ia^aj (33) 

< Re (Tr^~(L + h + ia)- 1 ^ < Re (T^(L + /i)" 1 ) < Re (Tr(L + /i)" 1 ) (34) 

|ImTV^(L + h + ia)~ l \ < - TV^(L + h)' 1 < - Tr"(L + h)~ l (35) 
here V\ is a generic positive measure subset of G. 

~cJ < /( a ) ( a (L + h- ia)- 1 ^ + h + ia^a) 1/2 (36) 



where /(a) < (a 2 (l max + h) 2 + a 4 ). 

These inequalities are easily obtained by applying the inequalities ( |32"U3"6| ) to the restriction to 
the Van Hove spheres Si r of L, a and h and taking the thermodynamical limit by letting r — > oo. 

□ 

Now we will prove the main theorem showing the existence of the bound for the magneti- 
zation. 
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Theorem 1 In an infinite and infrared finite graph G, Ve > 3/3 > such that V/3 > (3: 

lim M fc) > 1 - e (37) 
/i — >o 



We consider, in the thermodynamic limit, the expressions ( |29|) and (p0|) for M{(5, h) and for the 
global constraint, with the measure given by ([31]). At the end of this section we will prove the 
following lemma: 

Lemma 4 Consider a positive measure subset V\ ofV(G). If G is infrared finite, the quantity 



Xv x S'p >h (38) 



where S'g h {i) = jj^-Sp^ is bounded for all V\ when (3 — > oo only if the cti satisfy the condition: 



a 2 = 0. (39) 
This condition is true also in the limit h — > 0. □ 

We will call the set of values a» = {a\,a2, ■ ■ ■} satisfying (|39^ ) the stationary point of the 
statistical weight dfip^a). Now we separate the space of integration variables a, into two 
subsets Tpft(a) and its complement Fg^a) such that Fp^(a) is the region around the stationary 
point where the real part of the measure duplet) is positive. We define: 

e(/3,h) = \ l_ dRe(// i9> h(a))(I + h + ia)- 1 + - I cBm.(pp th (a))jL + h + ia)- 1 

z Jr p , h z J 

e'{P,h) = \ [ dRe(u8.h(a)) (%&(L + h + ia)- 1 + h 2 (L + h + ia)- 2 ] (40) 
z Jt^h VP / 

+~ J dlm(//g, ft (a)) + h + ia )^ + h 2 (L + h + ia)- 2 ^ 

Lemma 4 and the saddle-point theorem imply that: 

V e > 3/3' G R such that V /?>/?': lim \e{(3, h)\ < e/3 and lim |e'(/3, h)\ < e/3. (41) 

/i— »0 h-*0 

Now for M(P, h) we have: 

M((3, h) = € (p,h) + ~ [ dRe{^ h (a)){L + h + ia)- 1 (42) 

Z Jr^hia) 

= e{(3, h) + - [ dRe(nB,h( a )) Re (h{L + h + ia)- 1 ) (43) 

Z JTp, h {a) V ' 

The imaginary part of the second term in ( f42| ) does not provide any contribution to the integral 
because M(/3, h) is a real quantity. In analogous way for the constraint (|^) we get: 

1 = € 'm, h) + ~ [ dRe(fi h (a))Re (^Tr(L + h + ia)' 1 + h 2 (L + h + ia)- 2 ) (44) 
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Now in ( f43| ) we deal with a real quantity averaged with respect to a positive measure and we 
can use the inequalities of Lemma 3. In particular applying (|32]) in (|43|) we get: 



M(/3, fc) > e(J3, h) + ~ [ dRe(na h (a)) Re (/i 2 (L + /i + 



ia) 2 



With this inequality and the expression for the global constraint (E3j) we get: 

M{p, h)>l- e'((3, h) + e(/3, fc) - / da Re^ h (a))^Re (lF(L + h + ia)" 1 ) . (45) 

•/r ftfc (a) P v y 

Now we take the limit h — ► and we study the terms of this inequalities for large /3. For e(/3, /i) 
and e'((3,h) the behaviour is given by (41). For the integral term in (45) we can exploit the 
inequality (|34|) getting: 



< / da Re(/Wa))^Re (Tr(L + /i + ia)" 1 ) < ^Tr(L + h)' 1 f da Re(/Wa)). 

•/r flA (a) P v ' P -/r^ )fc (a) 

The saddle-point theorem implies that: 

3/3" 6 R such that V/3 > /?" : lim - / dRe(jiB h (a)) > 1/2. 

fe^O Z Jr Pth (a) 



(46) 



If for each e we choose 



we get: 



> (3" and > -2n lim Tr(L + hy 1 = (3"' (47) 

lim ^ / dRe(/^(a))Re (tF(L + h + fa)- 1 ) < e/3. (48) 

Finally fixed < e < 1, if we choose /3 = max(/?', /?", /?'"), we have that \/[3> (3 M(j3, h) satisfies 
condition (|37|). Note that from ( |47|) we have that 0" and then (3 are bounded only for infrared 
finite graphs. □ 



4.1 Inequalities 

Now we will prove the inequalities (p2|), p3]), (p4|), ( p5| ) and (p7| ) for finite graphs. On a finite 
graph the set of all real functions of the vertices (f>i has a vector space structure. Let us define 
{4>\ = (4>i, . . . , 4>n) and \(j>) = (0i, . . . , 0jv)*. The scalar product (0|^) is defined by: 

The matrices L, A, h and a are operators on this space. The matrix L is diagonalizable by a 
real transformation, and its eigenvalues I satisfy < I < l max - Let us prove the following: 

Lemma 5 On a finite graph, it is possible to introduce a real operator B defined by: 

B\L + h)B = I B l aB = c (50) 
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where c is a real diagonal operator and I is the identity operator. Furthermore B satisfies the 
properties: 

{L + h + ia)- 1 = Bil + ic^B 1 (L + h)- 1 = BB t (51) 
\\ Bt B\\ = \ ( 52 ) 

where 

..... (WBU) 
\\BrB\\ = sup- 



The existence of B is proved in |6j where it is also shown that is real. Properties (|5l|) can be 
immediately obtained by (|50|). Therefore it is easy to get the exact expression for B: 



B = TAT' (53) 

where T is the orthogonal transformation that diagonalize L; A is the matrix A^m = (1/V^fc + h)°~km, 
where Ik is the eigenvalue of L relative to the eigenvector k; finally T' is the orthogonal operator 
that diagonalize the symmetric matrix A^aTA. B is not an orthogonal transformation but we 
compute its norm: 



„ , „ WB'BUp) {4>\A Z \6) 1 

B'B = sup m IW = sup m -p = - 

This proves property ([52]). □ 



Proof of Lemma 1 Let be a generic vector of the space of the functions of vertices. We 
have: 

Re{<j)\h 2 {L + h + ia)-^) = Reh 2 {<f)\B(l + ic) _1 £*B(l + ic) -1 B*|0> 

= /i 2 (<£|B(Re(l + ic) -1 ).B*B(Re(l + ic)" 1 )^) 
-/t 2 (0|B(Im(l + ic) _1 )B*S(Im(l + ic) -1 )^^) 

< /i 2 (0|£(Re(l + ic)~ 1 )B t B(Rs{l + ic) _1 ).B*|^) 

< /i(0|5(Re(l + ic) -1 )(Re(l + ic)- l )B l \cp) 

< /i(0|5(l + c 2 )- 2 5'|0) 

< h^Bil + c 2 )-^ 1 ^) 



where we used properties (51) and (|52]). Then we have: 

Re(<f>\h(L + h + ia)- 1 ^) = h{cp\B{Re{\ + ic)- l )B l \cp) 

= hi^Bil + c 2 )-^^) 
< h^BB 1 ^) < 

so we get: 

Re((/)\h 2 (L + h + ia)~ 2 \(j)) < Re((f)\h(L + h + ia)- 1 ^) < (0|0). (54) 
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Now choosing in (54) (<fi\ to be the constant vector on = 1, we prove (p2|). 
To prove the inequality (23) we consider the following expression: 

< {(pKL + h-ia^iL + h + ia)- 1 ^) = {cf>\B(l - *c) -1 B t B(l + ic^B*]^) 

< h^Bil-ic^il + ic^B 1 ^) 
h 

< 1(01^(1 +c 2 r i B t \ ( f>) 

< ^Rei^Bil + ic^B 1 ^) 

< -Re(0|(X + h + ia)' 1 ^) . (55) 
h 

where we used © and (H). Putting <$> i = ati in ([55|), we get (|23f . 

To get inequalities ( p4] ) and (25) we introduce the base (/s|, k = 1, . . . , N in the vector space 
of the functions of the sites given by: 

f h=l ifi = * 

y in all others vertices. 

We have: 

(L + h + a)ll = (k\{L + h + a)- 1 ^) = (k\B{l + ic)- x B*\k) (56) 
For the real part of (|56|), with properties ( |5l|) we get: 

< Re(L + h + a)^l = (k\B(l + c^B^k) < (k\BB l \k) = {L + h)^ 
If we sum over all k E V\ we get: 

< ReTry A (L + h + a) -1 < Try A (L + /i)" 1 < Tr(L + hy l 
This proves (E5)- Now for the imaginary part of (|56|) we have: 



\lm{L+h+a)^\ = KA^l+c 2 )- 1 ^)! < (A^lcKl + c 2 )- 1 ^) < ^(k\BB t \k) = ±(L+h)£ 
If we sum over all k E Vx we get: 

|ImTry A (L + h + a)- 1 ] < £ |Im(L + h + a),", 1 ! 1 Try A (L + h)" 1 < ^Tr(L + h)~ l 



In this way we proved inequality (pq). 
For (^) we have: 

a 2 = (a|a) = (a|(L + h + ia)(L + h + ia) _1 |a) < 

i 

< (a\(L + h + ia)(L + h- ia)\a} 1/2 (a\(L + h- ia)' 1 ^ + h + ia)- l \a) 1/2 

( V /2 

< /(«) ^ Qii(L + h — ict) ik (L + h + za) fcl - a 7 - =0 

V iftj / 

where /(a) = (a|(L + /i) 2 + er^a) 1 / 2 < [E;(a 2 (^ mm + /i 2 + a^)] 1 / 2 and this completes the proof 
of Lemma 1. □ 
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4.2 The saddle-point condition 

Proof of Lemma 4- Let us compute the quantity (|38[). We have: 

^S^~ h = ~ Tt Vx (L + h + ia)- 1 + | ( M (F A ) - h\L + h + ia)-i X (Vx)(L + h + top 

where x(Vx)ij = Xv A <%■ The first term in this expression is always bounded, also when h — ► 
for infrared finite graphs; inequalities and (|25|) can be used for the proof. Then for (3 — > oo 
( |38| ) is bounded only when the condition: 



M(Vx) = /i 2 (L + /i + ia)" 1 *^)^ + /» + ia)- 1 (57) 
is satisfied. Let us consider the particular case where we put in fl57D V\ = G. We get 

1 = h 2 (L + h + ia)- 2 . (58) 



Let us show that (|58j ) is satisfied only when ( |39[) holds. First of all we have that (^) implies: 

Re (h{L + h + ia)- 1 } = 1. (59) 
Indeed with the inequality ( |32[ ) we get: 

1 = Re (/i 2 (L + /i + to)- 2 ) < Re (h{L + h + ia)- 1 ) < 1. 



Now using the Schwinger-Dyson identity (|3|) in ( |59|) we have: 



1 = Re (7t(L + /i + ia)" 1 J = 1 + Re ((L + h + to)- 1 (-to 

and then: 



= Re((L + /i + top(-to) 



Re ( + ^(-to)(L + /i + to)- x (-to) 
1- 



= Re + /t-HapoJ . (60) 

where in the first step we used again the Schwinger-Dyson identity. Now applying the inequality 
to ( pC| ) we get: 

< a{L + h — ia)- 1 (L + h + to) _1 a < Re ( — aj(Z + h + to) _1 0!^ 



a(L + /i - to)" 1 ^ + h + to)- 1 ^ = (61) 
Let us consider the inequality (j36[). We have: 

1/2 



n 



2 < /(a) (a(L + /i - ia)-!(L + /i + to)" 1 ^ = . 
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Notice that f(a) is bounded if a 4 is bounded. So we proved that the condition (39) for the 
variables a>i must be satisfied, if the quantity ( |38| ) is bounded when f3 — > oo. This proof also 
holds when h — > 0. Indeed in this case the stationary condition (j5q ) can be expressed by (|6l| ) 
and the inequality (^6|) holds. 

Now in order to complete our proof we must verify that also for a generic V\ (38) is bounded 
at the stationary point when j3 — > oo. In this way we will show that (|39] ) is not only a necessary 
condition but also a sufficient one. If we evaluate ( j38|) with the constraint (^9|) we get: 

/? 



2 -Trv A (L + h)~ L + | ( M (V A ) - h (L + h)-\{V x )(L + h)' 

-^(L + ^ + ^M^-xvT 

where we used again the Schwinger-Dyson identity. For infrared finite graphs the last expression 
is always bounded (also when h —* 0), see (pd|). This completes our proof. □. 
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